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Abstract 

We analyze a distinctive mechanism for inflation in which particle production slows down a scalar 
field on a steep potential, and show how it descends from angular moduli in string compactifications. 
The analysis of density perturbations - taking into account the integrated effect of the produced 

particles and their quantum fluctuations - requires somewhat new techniques that we develop. We 
then determine the conditions for this effect to produce sixty e-foldings of inflation with the correct 
amplitude of density perturbations at the Gaussian level, and show that these requirements can be 
straightforwardly satisfied. Finally, we estimate the amplitude of the non-Gaussianity in the power 
spectrum and find a significant equilateral contribution. 
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1 Introduction 



Inflation [1] is a very general framework for addressing the basic problems of primordial 
cosmology. It requires a source of stress-energy which generates an extended period of 
accelerated expansion. This can arise in many different ways even at the level of a single 
scalar infiaton, for which the space of inflationary models has been usefully organized by 
an effective field theory treatment [2]. The various mechanisms can be distinguished in 
many cases via their distinct predictions for the CMB power spectrum and for relics such 
as cosmic strings that may be produced after inflation. As well as being observationally 
accessible, inflationary theory is also sensitive to the ultraviolet completion of gravity, for 
which string theory is a promising candidate. 

Traditional slow roll inflation requires a flat potential, which can be obtained naturally 
using approximate shift symmetries or with modest fine-tuning. Inflation, however, does not 
require a flat potential. Rather, in general in single-field inflation [2, 3] the infiaton may 
self-interact in such a way as to slow itself down even on a steep potential as in e.g. [4-7]. It 
is interesting to examine such mechanisms further, to explore their novel dynamics and to 
better assess the level of fine-tuning required to obtain them from the point of view of both 
effective field theory and string theory. 

In this work, we analyze a simple mechanism for inflation in which the infiaton rolls 
slowly down a steep potential by dumping its kinetic energy into the production of other 
particles Xi (plus appropriate supersymmetric partners) to which it couples via interactions 
of the form 

l9'J2(<i^-<p^rx^. (1.1) 

i 

As rolls past each point 0i, the corresponding Xi particles become light and are produced 
with a number density that grows with increasing field velocity (p. As it dumps energy into 
the produced particles, slows down; meanwhile the produced particles dilute because of 
the Hubble expansion. With sufficiently closely spaced points Xi ^6 will see that this yields 
inflation even on a potential that is too steep for slow-roll inflation. This mechanism, trapped 
inflation, was originally suggested in [8] based on the preheating mechanism developed by 
[9] It can be usefully viewed [5] as a weak-coupling analogue of DBI inflation (or vice 
versa) in which the effects on 0's motion from the production of the x fields dominates over 
their loop corrections to its effective action. 



From the low energy point of view, although couplings of the general form (1.1) are 
generic, the prospect of many closely spaced such points (pi seems rather contrived. However, 
we will see that just this structure descends from string compactifications in a rather simple 



^ There are other interesting approaches using a gas of particles to slow the field evolution on a steep 
potential in order to inflate (see e.g. the recent review [10] and [11]) or to avoid the overshoot problem in 
small-field inflationary models (see e.g. [12-14]). The change in the CMB power spectrum from a single 
particle production event was also studied in [15]. 
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way. It arises in the same type of angular directions in field space that undergo monodromy 
from wrapped branes as studied recently in [16, 17]. 

In [16, 17], a single wrapped brane was considered. A scalar rolls down the potential 
over a large distance A0 ^ Mp corresponding to multiple circuits of an underlying circle 
around which the brane tension undergoes monodromy. In this super-Planckian regime, the 
potential satisfies slow roll conditions as in chaotic infiation [18] (though with a distinctive 
power law behavior depending on the example). In the same direction, at sti^-Planckian 
field values cf) < Mp, the potential is too steep for slow roll infiation. However, in variants 
of these setups, because of the underlying small circle, the system periodically develops new 
light degrees of freedom as the infiaton rolls down the steep part of the potential. 

The analysis of the perturbation spectrum in this class of models, including the integrated 
effects of the produced particles, has interesting novelties. The number of produced particles 
fiuctuates quantum mechanically, leading to a source term in the equation of motion for 
the perturbations of the infiaton. A constant solution to the homogeneous mode equation 
develops parametrically before the mode stretches to the Hubble horizon, as in previous 
examples of single field infiation in the presence of a low sound speed [3,6, 19]. Finally, as 
with other mechanisms such as [5,6] in which interactions slow the infiaton, a simple estimate 
reveals a correspondingly large non-Gaussian contribution to the perturbation spectrum in 
trapped infiation, which will be within the range tested by the upcoming Planck satellite [20] 
according to preliminary estimates for its capacity to detect or constrain the three-point 
amphtude f^^^^^^'^^K 

While this work was in completion we received the interesting work [21], which has some 
overlap with the present paper. 

2 Background Solution 

In this section we will find the background solutions and the conditions for trapped infiation, 
without making use of the perturbation spectrum. We will discuss perturbations in the 
next section. Getting the power spectrum to match observation will further constrain our 
parameters. 

The idea of trapped infiation is that particle production will slow the infiaton (0) enough 
to produce infiation on a potential which would be too steep for slow-roll infiation. For this 
to work, we will need a number of different fields to become massless at regular intervals 
along the direction. A Langrangian describing such a configuration can be written as 

C = ld,<pd^<P - V{<P) + \ Y.^d.Xid'^Xi -9\<P- <P^fXi) + ■■■ , (2.2) 

i 

where the . . . represent the supersymmetric completion of these terms, applicable in appro- 
priate cases. Softly broken supersymmetry helps to suppress Coleman- Weinberg corrections 
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to the effective action arising from the loops of hght x particles. As discussed in [5,8], 
at weak coupling particle production dominates over quantum corrections to the effective 
action for colliding locally maximally supersymmetric branes. Here the points 0i are the 
points where Xi become massless. For simplicity, we take these to be evenly spaced in with 
spacing A = — 0j. The coupling g may be small. If starts rolling down the potential 
V{(j)), whenever it crosses a point 0i, Xi particles are produced. The expectation value |^ of 
the number density of the Xi particles produced is given by [8,9] 

where ti is defined by 0(ti) = 0i and the powers a = accounts for the dilution of 

particles due to the expansion of the universe. The energy density of the x particles is then 
given by g\(f) — (t)i\ny^ following the particle production event, i.e. once the system has passed 
back into the adiabatic regime where ^ -C 1. Because m-^ = g\(j) — (f)i\, the x fields behave 
adiabatically when 

. (2.4) 

Making the replacement |0 — 0j| ~ (pdt, we can estimate the timescale on which the particle 
production happens: At ~ {g(p)~^^'^. Requiring this timescale to be short compared to 
Hubble At -C if-^ imphes 

<^g4) . (2.5) 



The parametric dependence of eq. (2.3) can be understood by noticing that the particle are 
effectively massless at production time, and are produced during a time of order At. This 
explains why n ~ 1/At^. On a longer timescale, Hubble dilution becomes important, and 
n oc a{ty^. 

The equations for motion for the homogeneous background solution (including the energy 
density in x particles) can be derived either from D^Tj^ = or by approximating XX with 
{xx) ill the equations of motion for as explained in [9]. The equation of motion is 

+ 3i70 + V'{<P) + jl^M^^))''"^^ = ' (2-6) 

where V = This sum over particle production events will be difficult to deal with, so 
we would like to replace it with an integral, giving us 

+ 3i70 + n0) + f ^M''))^'"^d'' = ^- (2-7) 



^For the purposes of calculating the homogeneous background inflationary solution, the expectation value 
of is all we will need. In calculating the perturbation spectrum in the next section, we will require its 
higher point correlation functions. 
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This is a good approximation to the sum when the variation of the integrand is small between 
production events. This is quantified by the two conditions ^ <^ 1 and ^ ^ 1. Because 
of the exponential suppression and the slow variation of the integrand, we can replace the 
integral by 



m)i"^d^ ^ ^.^mi ^ (2.8) 



roximation under the condition <^ 
that the particle production is the dominant mechanism for damping (|0| <^ 3i/|(/)| <^ V 



This is a reasonable approximation under the condition ^- If we assume slow roll and 



then we can solve (2.7) to get 

. (3i7A(27r)3r)i 
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(2.9) 



It is worth commenting on the limit if of the above expression. In this case goes 
to zero. This is due to the fact that in absence of dilution, the mass of the particles increases 
as (f) moves after the time of particle-production. Therefore, (j) loses energy even after the 
particles stop being produced. This explains why, in this H —>■ limit, the solution is 
different from the (p = const, that one would naively expect in the case of a linear potential. 
In the presence of a non-zero H, the growth in mass of the particles is compensated by their 
dilution, which allows for a steady solution (p ~ const, to exist. 



Given the solution for (p in eq. (2.9), we can find H and the slow roll parameters. The 
usual Friedmann equation is 

SM^pH' = p^ + p^ = + V{<P) + J29\<P- 0>x - yi<P) ■ (2-10) 

i 

We are assuming that the energy density is dominated by the potential energy. Using 
energy conservation p = —3H{p + p) we get 

QM^HH = -3H{<j)'^ + ^^|(/. - (j),\n^) , (2.11) 

i 

where we have used ~ 0. The generalized slow roll parameter e = — is then given by 

^_ 3(0^ + E.^?l0-0.K) ....^ 
^~ 2V ■ ^ ^ 

As expected e <^ 1 is the statement that the energy density is dominated by the potential. 

We would like to use e to constrain our parameters. We will assume that p^ ^ 0^ so our 
condition becomes 
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In order to constrain our parameters, we will make some estimates of this integral. Using 
— — ~ '^') (given |0| <^ \4'\H), we can do the integral to get 
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Using (2.9), V = 3MpH^ and dropping order one factors we get 

{2rr f'^—^ « 1 • (2-15) 
gMj,H^ 

We are now in a position to massage some of our previous inequalities to get conditions 
on individual parameters. Using (2.5), we can use our solution to get the inequality 

i/t < y'iAi(27r)^/^ (2.16) 

This provides a lower bound on A. The requirement that the particle production events 
were frequent also gave us the inequality A ^ \(j)\H^^. Using (2.9), this gives us 

At « ^i2nf/' . (2.17) 

These two inequalities imply 

gH^ < (27r)V . (2.18) 

We can also use our constraints on (p to get analogues of the slow-roll condition rj <^ 1. 
Recall that our solution required -C 3H(j) <^ V. Taking a derivative of (2.9) we get 

0=^(-ei/0 + 02^). (2.19) 

The first inequality, <C SHc/) is trivially satisfied for the first term, but the second gives us 
a new condition 

(2vr)6/5^-^ « 1 . (2.20) 
The second inequality, 3H(j) <^ V also gives a non-trivial condition 

(2vr)^/^^ « 1 . (2.21) 
gV 5 

There is another important requirement that we have ignored. Infiation is required to 
last long enough to give at least 60 e-folds. We will discuss this constraint in the context of 
an m^0^ model, after we discuss perturbations. 



■^In general, we will not keep track of all order one factors, in part because our analysis of the integro- 
differential equation governing (j) and its perturbations will not be exact. 
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3 Perturbations 



3.1 Gaussian Perturbations 

Determining the form of the curvature perturbation is a dehcate task. Since the trapping 
is intrinsically a multifield effect, we have not developed a Langrangian description of our 
effective equation of motion for cj) that one can consistently perturb. The strategy that we 
will use instead is to study the perturbations using the equations of motion directly. 

There are two approaches one could take. The first is to use constant 0, 'unitary', gauge 
and perturb in the metric. This would seem to have an obvious advantage, given that the 
particle production would happen everywhere at the same time in this slicing. Unfortu- 
nately, solving the many constraint equations for the metric perturbation is a complicated 
task. Similarly to what occurs in spontaneously broken gauge theories when one works in 
unitary gauge, this would also be the gauge where the main physical degrees of freedom are 
most obscure. As is usually the case in infiation [2], the matter scalar degree of freedom 
produces some scalar perturbations on the metric. These are not independent scalar degrees 
of freedom, but they are constrained variables. These perturbations of the metric are less 
important than the matter scalar excitations (the scalar field and the x particles in our 
case) for all the range of energies that we are interested in: from deep inside the horizon 
to freezeout. Thus, it is convenient to work in a gauge where the scalar field and the x 
particles appear explicitly, so that one can neglect the metric perturbations. This leads us 
to the second possible approach to study the perturbations, which will be the one we take 
here, where we work with constant curvature slices. In reality, since metric perturbations 
are less important, we will forget about them from the start, and we will work directly in an 
unperturbed quasi de Sitter universe. We will therefore perturb our equation of motion for 
0, taking into account the variance An(x,t) in the number density of x fields created. After 
horizon exit, at the time of reheating, these are converted to the curvature perturbation in 
the standard way [22]. 

The equation of motion for 0(x, t) takes the form 



^ + 3/70+^(0) + / + (3.22) 



We have assumed as before that we can make the sum of sets of produced x particles into 
an integral (an approximation to be checked below), and we have included their quantum 
fiuctuations in the last term. 



In the Gaussian approximation, the last term in (3.22) is equivalent to the variance in 
the number of produced x particles: 

9'Y.^x]-{x]m-h)^9^<x,t) . (3.23) 



7 



This behaves as a source term in the equation for the inflaton perturbations. This is some- 
what analogous to the equation for perturbations discussed in [23], and we can use some of 
the same techniques. We will now perturb the field around the background solution as 

0(x,t) = + ¥?(x,t) . (3.24) 



When expanding our equation of motion in if, we have to be careful to keep all the 
contributing terms. In particular, fluctuations of the inflaton change the time when particle 
production occurs at different spatial points. This manifests itself as a fluctuation of t', our 
variable of integration, when we use the continuum approximation to the sum over particle 
production events. We deflne t' by = (j){x,t'). Expanding in ip and t' — > t'^ + 5t\ we flnd 



5t' 



(3.25) 



We should think of the integral as being over tg. This implies that the upper limit of the 
integral is also subject to the perturbation. In particular, we can think of the integral as being 
over all time with a step function Q{t—t'Q — 5t'). This accounts for the fact that, on equal time 
slices, at different spatial locations, a different number of x flelds could have become massless 
and therefore been produced. In these regions, a different number of particles contribute to 
the sum, leading to a different region of integration. 

Putting these pieces all together we get the equation of motion for the fluctuation 



(p + —^ + 3Hip + V'{(f) + (p) + 



t-St' 



9' 



A(27r)3 

-9'J2[(^"-(x',m+v-<Pj)], , 



{<P{t' + 6t') + ^{f + 5t'))l^i^l±^dt' 



a{ty 



(3.26) 



where we have done a Fourier transform in the spatial direction with ^ being the physical 
momentum. For the Gaussian fluctuations, we will expand to linear order in ip. This gives 
an effective equation of motion 



a(t' 



./\3 



a{ty 



-dt' 



where we have deflned 
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A(27r)3 



-gAn{k,t) , 
(3.27) 

and used (3.23). One can check that V" <^ is the 



r^-like condition (2.20), so we will drop the V" term. We will see that this condition, not 



V" <^ H^, is sufficient to ensure that the spacetime is accelerating, and the modes freeze-out 
and produce curvature perturbations. This is very different than in the standard slow-roll 
case. 

There are two types of contributions to the power spectrum - those sourced by An, and 
those which would arise in the absence of the source. We will flnd that the former dominates. 
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To begin, in order to analyze both these contributions, we require the homogeneous mode 
solutions to the above integro-differential equation. This will allow us to construct the 
Green's function required to determine the sourced perturbations. 

To get some intuition for the behavior of the homogenous solutions, we will start by 
solving the equation for constant p = ka~^. This is a good approximation when <^ 1 
which holds until p c:^ H. We will also approximate H and as constant, which holds to 
leading order in our generalized slow roll parameters. There are three epochs of interest 
depending on the ratios p/yh and p/H: 

(I) p ^ m: The modes are approximately Minkowskian, with both Hubble friction and 
particle production effects negligible in their equations of motion; we start with the pure 
positive frequency modes corresponding to the standard Bunch-Davies vacuum. 



(II) H <^ p <^ m: In this regime, a constant solution to (3.27) appears. The mode solu- 
tions from region I, evolved into region II, develop a term which is approximately constant. 
This contribution begins with a very small amplitude (which will be determined in our ex- 
act solution below) but ultimately dominates over the other terms which become damped 
exponentially in Ht. 

(III) p < H: In this regime, the curvature perturbation ^ = ^Lf becomes constant, lying 
outside the Hubble horizon. 

In particular, we will find that the modes actually freeze-out well before reaching the Hubble 
horizon. This is somewhat analogous to the freeze-out of modes at the sound horizon Cg/H <^ 
1/H in general single field models of inflation [4-6, 19]. 



Now let us derive these features from a more detailed analysis of (3.27). For constant p 



rh and H, we can find exact solutions to (3.27) using the ansatz (p{k,t) oc e"*. We can solve 



the equation trivially because in our WKB regime of constant p, all terms are proportional 
to e"* with constant coefficients depending on H and a. In particular, using the ansatz and 



doing the integrals we find that (3.27) reduces to 

5 



(a^ + 3Ha + p^ + rh^ + = . (3.28) 

This equation gives the mode solutions when a ^ —3H. Multiplying through by 3H + a, 
we get the cubic equation 

+ 6Ha^ + {9H^ + + m^)a + 3Hp^ = , (3.29) 

where we have defined rh^ = l^rh?. It should be clear from this equation that behavior of the 
perturbations will only be different from the usual case if iv? ^ H^. In this model, this is 
always the C3iSG, clS this condition is equivalent to the slow roll condition 3H(j) -C V . 



There are three analytic solutions to (3.29 ) since it is a cubic. To understand the behavior 



of the solution and impose boundary conditions, it will be useful to expand these solutions 
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perturbatively in the different regimes discussed above. When H"^ <^ P^, we can expand the 
modes around H = 0, giving 

^- 3H2rh^+p'^ p2 

a± = ±2 VP - — — 5- ) «3 = -3f/ — - — ^ . (3.30) 

When ^ we can match onto the solutions in the Bunch-Davies vacuum. Specifically, 



we should use the mode with a normalization of l/i/Sp. Using (3.30), taking into account 
that e"+* dies exponentially like l/a^^^, we see that this corresponds to a Minkowskian mode 
solution of the standard normalized form 



k 



u^it) = -^e-^^ . (3.31) 
ay K 

When p drops below m we need to match onto the modes in the rh ^ p regime. Notice that 
in the limit, the mode decays very slowly compared to a±. In essence, these modes have 
frozen out at the scale k = ma. 

One might have worried that when V" > H'^, the fluctuations of would be massive and 
would not produce curvature perturbations. Like in small speed of sound models, we find that 
the mass can be much larger than H and still contribute to the power spectrum. Replacing 



p^ ^ p^ + V" in (|3.30[), we find that ~ -?>HV" /ni^. Therefore, as long as V" < 



there is still a nearly constant mode that will be converted to curvature perturbations. This 



condition is equivalent to (2.20) and is always satisfied in these models. 

When matching the modes at p ~ m, it is clear that the leading terms in a± are smooth 
at the cross-over. The real part of a±, however, transitions from — 3if/2 to — 3if in the 
crossover between regions I and II. This behavior is distinct from what would arise for a free 
scalar field in de Sitter space, and the matching between the two solutions will introduce 
new effects suppressed at small 0{H/m). In order to determine the relative amplitudes of 
the modes, we cannot simply match the two regimes using continuity at p ~ m. Such a 
matching calculation assumes that crossover is rapid, but the wavelength ifr^ of the modes 
at the crossover is much smaller than the time period during which the crossover takes 
place. Therefore, in order to calculate this sub-leading contributions we will need more than 
the WKB mode solutions. 

Let us therefore move on to discuss the exact solution to the homogeneous linearized 
equation for the perturbations. It proves to be convenient to transform the equation to 
conformal time r = —1/aH, with late times corresponding to r — ^ 0. Denoting the derivative 
with respect to r by ', we have 

,, 2 , , r, m^r r dr' f b , , \ qAn(k, r) 

- -y + k'v + + ^ / 37 o^V + 3^. = -'—h^ . (3.32) 



Let us comment on the structure of the source on the right hand side of equation (3.32) 



Since the particle creation happens on very short time scales, we can concentrate on the 
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Minkowski limit. In this case, the squeezed state describing the created x particles in the 
case of homogeneous motion takes the form 



(3{kp)a}-. a) 

tip 

2a* {kp 



1^) = A/'exp 



E 



|0) 



(3.33) 



where a, (3 are Bogoliubov coefficients satisfying |ap — = 1 and A/" is a normalization 
factor. Here kp represent the physical momenta, given by kp = k/a{t), where k is the standard 
comoving wavenumber. From t his, o ne computes the expectation value of the number density 
J cPkplPj: \/{2nY given in (2.3), using the standard result (reviewed in [8]) that 



(3.34) 



Similarly to the case of the computation of the expectation value of rii, where i represents 
the Xi particle species, it is quite straightforward to see that 



{An^{k,t)Anj{k',t')) 
{2rrf6^''\k + k')6i 



(3.35) 



3/2 



a(t)3/2a(t')3/2 



aitf 



a{t') 



Here the 0(t — tj) function (and analogously 0(t' — tj)) represents the fact that, for the 
population i, particle production is irrelevant before the particles become massless. This 
is only an approximate expression, which is parametrically correct but that we expect will 
receive order one corrections in a full calculation. The purpose of this first paper on this class 
of models is to understand the main features of the predictions, and therefore we consider 
this level of accuracy enough for the present. By using the definition 



n(A;, t) = ni{k, t) 



(3.36) 



we obtain: 



{An{k,t)An{k',t')) ~ ^(Ani(A;, t)An^-(A:', t')) 



{27if6'^^\k + k')J2 

i 

We can substitute as usual 



3/2 



a(t)3/2a(t')3/2 



e(t - 



a{ti 



aity 



a{ti 



a{t') 



E 



dti 



A 



(3.37) 



(3.38) 
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to find: 



{An{k,t)An{k',t')) (3.39) 
It is straightforward to see that the integral gives: 

where teariy.Uate are the smaller and greater of t,t' . Here A^hits ~ (f>/{HA) is the number of 
particle production events contributing; because of Hubble dilution, this is limited to events 
taking place within a Hubble time. 

Later in the section, we will see that the An fiuctuations source the infiaton perturbation 
through the integral in cosmic time of a Green's function whose width in time is of order 
H^^. This means that the infiaton perturbations will be sensitive only to the integral in time 
of the correlation function of An, and therefore we can approximate the time dependence of 
the above equation with a 5— function to obtain: 

{An{k,t)An{k',t')) ~ i^^{2n)'6^'\k + k')N^,,,H-'6{t - t') (3.41) 

We stress that this expression would receive order one corrections in a more exact calculation, 
but we expect it to capture the correct parametric dependence of the result. 

Finally we note that this expression can be obtained more directly in the case where 
there is a single production event per Hubble time (and correspondingly Nhus species in 
this time). Then, the particles from the jth event have diluted significantly before the next 
occurs, and the time dependence of the correlation function can be modeled approximately 
using J = tjH by NhusSij = Nhus6iH{t - t')) = NhitsH-^5{t - t'). 

It is convenient to rewrite (3.32) in differential form by acting on it with t^^, giving 

/// 4 , 3 ,, , 9 , iv? , k"^ / : N 3 /r- — — d Ah , , , 

V + -^V - -V + fcV + JT^^V ~-V = -9i9(l))-WNhUsH~W- , 3.42 

where An is defined to have unit variance: 

{Ah{7], k)Ah{7]', k')) = {2nf6^^\k + k')6{r] - r]'). (3.43) 
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In this form, the general homogeneous mode solutions (fhom can be written in terms of 
hypergeometric functions, expandable in terms of Bessel functions. We find (using Mathe- 
matica): 

/ 1 im 1 im 1 1,22 

ZH — irh ZH — irh l 

+ C2X2 T^k^^iF2\l 



+ C3 X 2 H k H -^Fr 



im 


5 

2 " 


im 
2H 


1 - 


im 


4 


+ 1; 


im 
2H 


5 

+ 2' 


im 
~H 


+ 1; 


--A;V 
4 



if 



2if 

= 5^c,/,(r). (3.44) 



i=l 



The function /i goes to 1 as r 0, and it represents the late time constant mode; the other 
solutions /2,3(t) decrease to zero as r — > 0. Imposing that this match the Bunch-Davies 
vacuum solution at early times yields three conditions on the three constants Ci,C2 and C3. 
We find that for large m/H 

Ci cx e 2ff . (3.45) 
This leads to a tiny contribution to the power spectrum from homogeneous modes: 

(g^ + m^)^ech^(fg) 1 

Because of this exponential suppression, the homogeneous contribution will prove to be 
highly subdominant to the sourced contribution. 



To calculate the perturbations generated by the source (3.23), we must determine the 
Green's function for the differential equation (3.42). We can define Gfc(r, r') as the solution 

to 

+ -,G', - ^-Gl + eC, + -^G; - -G, = 6{r - r') , (3.47) 



r 



with the boundary conditions across t = t' given by 

AGl = 1 , AGk = AG'k = . (3.48) 
It is useful to change variable from t to f = kr, and solve the simpler equation 
•~ 4 ~ 3 " - m^ ~ I ~ 

G + —G-^G + G + -—G-^G = 6{T-f'), (3.49) 
with the boundary conditions across r = t' given by 

AG = 1 , AG = AG = . (3.50) 
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Here a dot stands for a derivative with respect to r. Notice that in this way all the dependence 
on k is implicit in the definition of f, and we have the simple relation: 



^G{kr, kr') 
k"^ 



(3.51) 



This change of variables will allow us to see analytically that the power spectrum is scale 
invariant. The sourced perturbation is given by 



3 / 

Vk{.T-) = g{g(f))i ^/ NhitsH- 



T' ■ ■ dr 

and the power spectrum at late times (r = 0) is given by 



C 



g"^'/'NutsH-' / dr' 



k^ 



X 



7/2 



NhitsH" 



dr' 



1 d 
r' dr' 

1 d 
f' df 



7 (r'G.ir 
f'G{f 



0,r')) 



(3.52) 



(3.53) 



0,f' 



This shows that the power spectrum is scale invariant. In order to determine its amplitude, 
we need to perform the integral above in (3.53), where we see that the power spectrum is 
determined by an 'effective' Green's function 



1 d 



f'G{f,f': 



3 

E 

i=l 



(3.54) 



with only the fi{T)'ji{f') term surviving as f 0. Though we have an analytic expression 
for the 7j's |^ we are unfortunately unable to perform the integral analytically. However, 
we can notice that the function 7i(t), whose only parametric dependence is on m/H, has a 
peak at the point f* ~ —m/H (corresponding to a physical momentum p = k/a ^ m), with 
amplitude 7i(r*) ~ H'^f^/m^ and width This allows us to estimate the integral (3.53) 
and to obtain the power spectrum: 



Pr 



X 



10 



-9 



1 



(3.55) 



Am k^ k^ 
This expression can be verified numerically. 

Finally, we should ensure that our integral approximation was valid in this context. When 
p is large, it is clear that the variation of if is large compared to the spacing between particle 
production events. However, the contribution from the integral only becomes important 
when the frequency of the modes is m. Therefore, the integral is a good approximation 
when mA0~^ <C 1. This condition becomes, using our background solution n2.9\, 



— \ - 

g2/\l 



< 1 . 



(3.56) 



*That we do not reproduce here for the sake of brevity. 

^Notice that, as anticipated, in cosmic time, this width corresponds to a time interval of order H~^. 
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This is a stronger version of the constraints sketched after eq. (2.7). 

We derived the integral term assuming that the particle production at each point is the 
same as for a homogeneous field, which is a valid assumption when the modes of interest 
obey -C gcf). Since the integral term becomes important at the freeze-out scale ~ fh?, 
we must have gcj) ^ fh?, which gives 

(27r)i2/5A4/5 



^iJl/5\/'l/5 



> 1 . (3.57) 



This constraint is similar to (and stronger than) eq. (2.5), but the origins of the two con- 
straints are different. In the next section, we will look at how all these constraints fit together 
in a model with V = ^rri^cp'^. 

Before moving on, let us comment on the role of the x fields in the perturbation spectrum. 
Our model is not a single field model, given that we require many x fields in order to slow 
the inflaton. As such, one might wonder if these extra fields may contribute to the density 
fluctuations. This is not the case because their mass grows to be large before any modes could 
freeze out. Specifically, the effective mass of a x field is given by ?^e// — d'^l'Pi't) ~ A 
Hubble time after the field becomes massless, the effective mass is given by '^e// — g'^<P'^H~'^. 
One can check that '^e// 3> is equivalent to our constraint (2.5). As a result, the x 
fields are massive compared to the Hubble scale and do not contribute to the curvature 
perturbation |^ 

However, given the crucial role of the x's in both the background solution and the gen- 
eration of perturbations, one must ensure that interactions do not cause them to decay. By 



^It is interesting to consider the fate of these heavy particles. In some regions of our parameter space, 
they are always lighter than Mp: g{(j) start — 4>end) ^ Mp where (j) start and (jiend refer to the start and end of 
inflation. If in other regions they become heavy, they may decay (certainly Planck mass black holes decay 
rapidly to lighter species of particles) . 
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construction, the two-body decay XX ~^ is always present. We can ensure that none of 
our results are affected by this process by requiring that dilution of particles due to expan- 
sion is the primary cause of decreasing number density. This is expressed by the constraint 

•—1/2 

n{av) ^ H. Assuming cr oc g'^{4>{t — U))''^ and v gcj) {t — we get the condition 
{t — ^ g^{g4'H)~^ . Evaluating this expression at the moment the x fields are created, 

{tc — ti)^^ = \J~g^ leads to the constraint 

H->g^\JV^ (3.58) 

We will impose this constraint on our parameters although it is possible our results would 
not be significantly affected even in regions where it is violated. The mechanism itself can 
tolerate some production as long as the energy density from the decaying x's does not 
interfere with the perturbations. 



Let us also compare our result for the scalar power (3.55) with the curvature perturbation 
one obtains from the fluctuations in x energy density coming from the variance in x particle 
number on the right hand side of Einstein's equation. We can estimate this contribution as 

Ml^Cgravity ~ m^An . (3.59) 
Here (gravity is not the curvature perturbation ( but comes from the gofj, components of the 



metric. The expression (3.59) arises from the Hamiltonian constraint. This contribution is 
not directly contributing to a measurable power spectrum, but we would like to ensure that 
the curvature it induces during inflation is negligible. 

By going into Fourier space, and using the fact that the fluctuations are evaluated when 
fc/a ~ if we obtain, after using eq. (3.40): 



{CgravUymgravUyik')) ~ ^ X gV^ J^Nf,UsS^''\k + k') . (3.60) 

Notice that ~ gcpH^^ in this estimate. By comparing with the contribution we have just 
computed, ~ ^ip, we obtain 

> l2 / ~ \ 1/2 

^gravity _ /r { ^ \ 



This ratio (f^ /V is the usual slow roll parameter, which is much smaller than e in this model. 
The ratio fh/H has to be smaller than ~ 10 because of the constraint coming from non- 
Gaussianities (see next section). For the specific case we will study in the next section, where 



V{(j)) = m /2, the above expression is also equivalent to eq. (2.21) times an additional 



suppression from the number of e-foldings, and therefore it is always satisfied in that model. 



^This is due to the fact that the fluctuations An average quickly to zero on scales longer than H ^, and 
therefore the induced metric perturbations become constant after having redshifted up to the scale H. 
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3.2 Non-Gaussian Perturbations 



The size and shape of the non-Gaussian contribution to the perturbations are particularly 
important for distinguishing between different models of inflation [24]. Since our interactions 
slow the inflaton on a potential which would otherwise be too steep for inflation, we should 
expect a substantial non-Gaussian correction to the power spectrum as in [6]. A detailed 
prediction for the bispectrum requires the calculation of the three-point correlation function 
of the curvature perturbation, as first completed for single-field slow roll inflation in [25,26]. 



Following [23], we can expand the equation of motion (3.26) for the perturbation into 
first order, second order, and higher order pieces: 

= (^(1) + y,(2) + _ _ _ _ (3_62) 
It is again useful to translate the expanded equation of motion into conformal time, and derive 



its differential form (as done for the linearized equation in (3.42)). Then we can obtain the 
second order perturbation ^9*^^^ by integrating against the Green's function (jfe(r, r') the terms 
in the expanded equation of motion which are second order in Lp^^^ and An. By looking at 



eq. (3.26 ), one sees that one contribution comes from the expansion of the term proportional 
to in the equation of motion, giving a contribution to the second order perturbation of 
order: 

^'ik^) - I dr'G,ir,r')^{^'iT')^'iT')), . (3.63) 

Another contribution comes from taking into account the time delay of the perturbation 
inside the integral, giving rise to a term of the form: 

^S.(r) - J dr'G,{r,r')^{v"{T'MT')), . (3.64) 

Yet another contribution, of order Lp^^^An, comes from expanding the {g(j)Y^^ coefficient in 
the source term, giving a contribution 

(2) , ^ 7/4 /77 FTTI [ . ,n I '\ ' ^ (viW(r')An(r'))fc , . 

V'UnM - 9 ^ <P WNhitsH ^ / dTGk{T,T)T — . (3.65) 

There are additional terms coming from the expansion of m^, but it is easy to see that they 
give subleading contributions. Also the contribution from the non-Gaussian statistics of An, 
which in the absence of interactions can still come from x-particle shot noise, is expected 
to be negligible if the number of particles is large enough. This is in fact always the case. 
Estimating the size of the non-gaussianity of An by {fh^ /{n^NfiUs) Y^"^, where we have used 
that m is the typical scale at which the Green's functions peak, it is easy to see that in our 



model, by using eq. (3.55), this ratio is smaller than 10 °, corresponding approximately to 



a negligibly small /nl ~ 0.1 . 
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The three point function of our perturbations is of the form 



,(1) 



(3.66) 



and we are interested in this amphtude at late times, r 0. We can estimate this using the 
same method we used for the Gaussian power spectrum, and let us start with the term in 



eq. (3.63). The perturbations on the right hand side of (3.63) can come from any of the three 



modes (3.44), not only from the constant mode /i . This is so because the perturbations in 



Lf^^^ that source the second order Lp^"^^ in eq. (3.63) can be evaluated when still well inside 



the horizon when neither of the three modes has yet decaied. One of the leading effects we 
find comes from the /2 and fs modes giving a contribution to the three point function of 
curvature perturbations of order 

{2nf6^'\J2h) (I)'/ drCMr)^ [g{g^fW NhusH-^)' (3.67) 
X dT'gk^{0,r')g2.^{r,T') dr" gks{0,T")g'.^(T,T") + sjmm. . (3.68) 



If we pass to the Green's functions G, g defined as in the former section, we find: 
'hY f df G(0,f) fcim^ 



(27r)35(3)(^fc, 



i=l 



h k( fH(j) 



^(g{g^)i^NHUsH-^ 



(3.69) 



"'-''-df'mr')~92{rk2lh,f') f'^"^ df" ~g{0,f") -g'sjfh/Kf") ,^ ^ 

k2 I ; K3 + symm. 



k?, fca 



k. 



3 3 



i=l 



kf 



X2X3 



~ 2 

m 



/TX2 fTXi 
dfg{0, f)g^{fx2, f) j df"g{0, f")g^{fx3, f") + symm. . 



where we have defined X2 = k2/ ki and X3 = ks/ki. The former expression is of the form 

1 



l^jj [gig^f^^fhUF^) j;^g{x2,Xs,m/H) + sjmm. (3.70) 

3 

= (27r)35(3)(^^^)i7^(fc^^fc2,A;3) 

i=l 



The factor of which characterizes the dependence on the global scale of the momenta, 
tells us that the signal is scale invariant [24]. 



similar term with a pair of /2 or a pair of /a modes will change the final result by no more than an 
0(1) factor. We study the 72/3 term above because the cancellation between phases is particularly simple. 
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As we discussed above, the Green's functions G{0, r) or g{0, r) are peaked at r^, ~ 
—m/{Hk), and the product of integrals forming the r integrand also exhibit a peak at 
this value. We can estimate the size of F{ki, ^2, ^3) using knowledge of the peak at r* and 
series expansion of the Green's functions around r = 0. We find that the Green's functions 
g{0, f) and 5'2,3(t, r') can be expanded as a series of the form ^ a„(f/f^,)" with order one 
coefficients a^. Physically, we beheve this occurs because the only features of these functions 
occur near f ~ f,,,. Therefore, we only expect any non-trivial behavior when f ~ f^,. Using 
the above Taylor expansion at r ~ r* is likely inaccurate, but we think it should be reliable 
for order of magnitude estimates |^ 



In the language of eq. (3.54), the leading terms in the expansion for /2,3(t^) are of order 
jL3jL±i!g and those for 72,3(t) are of order ^^f^^^ . The oscillations contribute a suppression 
factor ^ to the integrals and an ^ enhancement to the derivatives. Altogether this gives us 
an estimate, which we checked against a numerical integration, of order 

3 tt2 7i\r2 

(c^.c,.Cg3) - (2-)^5(3)(5^ ^')^Mzr^ • (3.71) 

i=i ^ 

Here we were not careful with the momenta dependence, and the factor k denotes only the 
typical size of the wavenumber. A more careful numerical analysis for the shape function 
X2x\F{l,X2iX^) / F{1,1,1) as defined in [24] shows that most of the signal is concentrated 
on equilateral configurations. The equilateral shape can be understood to be a result of 
the Green's functions being peaked at a scale t^: we get a large contribution when all the 
momenta are equal and all the Green's functions can be evaluated at their peak value. More 



in detail, by looking at eq. (3.69), one can notice that in order for the integrals in f and f 
to include in their domain the peaks of ^(0,f') and of ^(0,f") by the time (7(0, f) reaches 
its peak at f ~ f*, we need to have 0:2, X3 < 1. However, in the limit f — > 0, we have 
approximately ^23(^2,3,^') oc ^^2:2 3, which suppresses this contribution to the shape by 
x\x\ and forces the dominant contribution to come from the case where X2-, x^ are as large as 
possible compatibly with the former constraint. We obtain that the integrals are peaked for 
0:2, X3 ~ 1, on equilateral configurations p°| The suppression of ^2,3('^5^') small f comes 
from the fact that the oscillating modes decay at late time. Notwithstanding the fact that 
the leading mechanism for generating non-Gaussianities is intrinsically a multifield effect, 
we conclude that the signal on squeezed configurations is not large, as is always the case in 
single field infiation [25-29]. 

(1') 

A similar analysis shows that the contribution due to '^\ si{t) is parametrically the same 
as the one of v^i.^^(T), while the one from ^^k\n{'^) suppressed by a factor of H/fh. The 



^AU the results using series expansions have been checked against numerical integrations and provide 
reliable estimates. 

^°There is some support also on flattened triangles, but the numerical study plotted in Fig. [2] shows that 
this does not dominate over the equilateral shape. 
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0.0 



Figure 2: A numerical study of the shape ^i^sP^O-^r^^^^s) choice of parameters ^ = 10, 

plotted in the region < X2 < 1; 1 — X2 < X3 < X2. The peak in the equilateral limit is 
clearly visible. 

remaining terms that we did not show are subleading as well. Since we are not careful with 
order one coefficients, there is no need to perform the calculation for ^plsti'^)^ since we do 
not expect cancellations or the shape to be peaked in the squeezed limit. 

Summarizing, following the standard definition, we can estimate the size /nl on equilat- 
eral triangles (with \ki\ = k) to be of order 

/■equilateral \^fci ^fe / /o ^r,\ 

/nl ~ -J^;^ ~ Jp ■ (3.72) 
where the primes indicate that we dropped the delta functions of momenta. 



4 The case V{(f)) = Im^cj)^ 

Let us now check the conditions for a viable model of trapped infiation, including the 
background solution and Gaussian perturbations. We will take a model with potential 
V{(f)) = |m^0^ for simplicity; other cases of interest include more general power law po- 
tentials Va{4>) = Given the number of e-foldings, the Gaussian power spectrum 



(3.55), and our solution (2.9), we can solve for two of the parameters and then express the 
various inequalities prescribed in §2 in terms of fewer model parameters. From this, we 
obtain the following relations. 

The number of e-foldings is 

fH,, 10-6 f MpY^'' f (P 10-6 



20 



where the slow roll condition (2.15) forced the total field range — 



= d) and we used 



(3.55) to eliminate A. Using this, (3.55), and our solution (2.9) to write the self-consistency 



conditions in terms of m/Mp, the constraints (2.16), (2.17), (2.21), (3.56) respectively give 
four conditions E3 



m 
Wp 



< 1 



< 1 



(27r)3l0 V/'iV3/2 f 



m 



m 



\Mp 

1/2 



3/2 



< 1 



< 1. 



(4.74) 



The constraint (3.57) and (3.58) are in the other direction: 



m 
Wp 



>1, 



m 
Wp 



> 1 . 



(4.75) 



Together, these constraints define a viable window in the space defined by the two free 



parameters [g, ■^), which is plotted in Fig. 



Maximizing the value of the field range 



Mp 



over this window, the field range is constrained to lie no more than an order of magnitude 
above the Planck scale. The mechanism therefore can operate below the scale typically 
needed for standard slow-roll infiation - in fact the constraints described above allow field 
values far below Mp, although we will see in the next section that experimental constraints 
on the size of the non-Gaussianity in the power spectrum prevent us from going far below 
the Planck scale in this model. 



5 Observational predictions 

In this section we will outline the predictions for the CMB derived from our infiationary 
mechanism. 

5.1 and r 



Because H and other background parameters change slowly during infiation, our power 
spectrum is approximately scale invariant. Its tilt is given by 



rio — 1 



dink 



Hdt 



(5.76) 



From (3.55) this becomes 



1 + 



H 



0.7 1 - 
1 — - ^ 0.99 



(5.77) 



^^The constraints (2.151, (2.201 give ^ 1, and so are trivially satisfied 
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Figure 3: Top: The allowed parameter window for the m (p model. The red zones are 



forbidden by eq. (4.74); the blue zones are forbidden by eq. (4.75) and by the constraint 



(5.82) on the size of the non-Gaussianities, to be discussed in the next section. The dashed 
line indicates the range of parameters for which 0/Mp ~ 1, with super-Planckian field 
ranges above and sub-Planckian ranges below. Bottom: Same plot as above for a model 
with potential equal to We do not give explicitly the constraints in the paper as they 
are very similar to the ones for the m^(j)'^ model. 
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where in the next to last passage we used eqs. ( 4.73[ ) and (3.53), and in the last passage we 
have used a typical number of e-foldings A^^e ~ 55 The parameter Q we have introduced 
here represents the ratio 4>cnd/<pNe, between the value of at Ne efoldings to the end of 
inflation and the value at the end of inflation. This parameter was not introduced in the 
former estimates because it does not affect them significantly. The condition e -C 1 requires 
that Q~^^^ — 1 -C Ne, so we can safely take Q ~ 1/4. The tilt is red, falling quite close to the 
statistically preferred region of the WMAP 5-year data [30] . However, it is worth mentioning 
that the tilt may not be a sharp prediction of this class of of models, but may be tunable in 
general. It can depend not only on the potential and the field range, but also on other details 
such as variation in the spacing between the particle production events, and in the mass and 
species numbers of of the particles. In order to compute the tilt very precisely, it would 
also be important to systematically check the contributions of higher dimension operators. 
These are limited by symmetries in our string-theoretic backgrounds, but we have not done 
a complete analysis of their leading effects. 

The power in gravity waves is as usual Ptensor ~ (^)^; leading to a tensor to scalar 
ratio of: 

ArfiH o 0-7 , f m \ 

g^NllQ^\2Tif — . (5.78) 



in the |m^(/)^ model. Maximizing this quantity over the allowed range of {g, j^) from the 
previous section, we find r <^ 10~^ for this potential. 



5.2 Non-Gaussianity 



Current constraints from data [30,31] bound j^'v^'^^*'^''^^ go that, using eq. (3.72), we have 



m 
H 



< 10 . 



(5.79) 



Note that since we have not been keeping track of 0(1) factors, there is a possibility that 
these may shift this constraint slightly in either direction. 

Plugging this into our solution, this is equivalent to 

V 



(27r)3^/73 



< 10^2 . 



For m this corresponds to 



1 [MpV [M, 



{2'K)'^g \ 4> J \m 



^ < 10 
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and imposing (4.73) we obtain 



m 10-1° 



(5.80) 



(5.81) 



(5.82) 



i^For potentials of the form - 1 = -(2 + 7a)(l - Q"^ /'^-°'/^°) /{2N^{U - a)) 
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This goes in the opposite direction from the previous conditions (4.74) (except for (4.75)), 
but leaves a wide window of viabihty. As seen in Fig. [3} it is the constraint on the non- 
Gaussianity that restricts the field range from going far below the Planck scale. This is to be 
expected - as in [5], [6], as the potential grows steeper a stronger interaction will be needed 
to slow the inflaton, and a larger contribution to the non-Gaussianity will be produced. 



6 Trapped Inflation from String Theory 

"Meetings are a great trap..." -John Kenneth Galbraith 

Because inflation is sensitive to Planck-suppressed operators in the effective field theory, 
it is generally of interest to model it in a UV complete theory of gravity. String theory, as 
a candidate UV completion of gravity, is a standard framework in which to develop such 
constructions. The present work was motivated in part by the top-down appearance of 
the structure required for trapped inflation. In this section, we will explain this structure 
and analyze the conditions for realizing trapped inflation consistently with moduli stabiliza- 
tion in appropriate examples. These realizations use the same structures recently used for 
monodromy-driven large field inflation [17,32], but now in a < Mp range of field. Because 
the relevant setups were described in detail in these works, our discussion here will be some- 
what more telescopic; the reader may therefore find it easiest to refer back to the relevant 
portions of [17,32]. 

To begin, consider wrapped D4-branes in type IIA string compactifications on nilman- 
ifolds, as in [16,32]. The simplest example of a Nil manifold suffices to exhibit our basic 
mechanism for closely spaced particle production events, though we will see that trapped 
inflation in this specific example would introduce too large a back reaction on the internal 
geometry. We will therefore ultimately be led to construct it in string theory by using axion 
moduli in warped Calabi-Yau compactifications of the kind analyzed recently in [17]. Particle 
production in these models was also considered in [33] where it was used for reheating. 

A nil 3-manifold is obtained by compactifying the nil geometry 

ds%ii = -^d'^i + PL^dul + Ll idx + — [uidu2 — U2dui]] 

= -^dul + (3Lldul + Ll{dx' + Muidu2f , (6.83) 

P 

(where x' = x — ^UiU2) by a discrete subgroup of the isometry group 

tx- {x,Ui,U2) {x + l,Ui,U2), 

. . . M . 

: (x,Mi,M2) [X - —U2,Ui + 1,U2) , 

M 

■ {x,Ui,U2) {x + —Ui,Ui,U2 + 1) . (6.84) 
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This manifold can be described as follows. For each ui, there is a torus in the U2 and 
x' = X — ^uiU2 directions. Moving along the ui direction, the complex structure r of this 
torus goes from r — r + M as Mi — Mi + 1. The projection by tui identifies these equivalent 
tori 
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At all values Ui = j/M for integer j, the two-torus in the U2 — x' directions is equivalent 
to a rectangular torus 

dslect = Lldyl + (3Lldyl , (yi, 7/2) = {Vi + ^1, 1/2 + ^2) , (6.85) 

(since r — r + 1 as j — j + 1). These coordinates yi and 1/2 are related to x' and U2 by 
an SL{2,Z) transformation. The 1-cycle traced out by U2 = A, A G (0,1) becomes a cycle 
(2/1, 2/2) = (MA, A) as Ml ui + 1. 

Consider first, as in [16], a D4-brane wrapped on this cycle. Near mi = 0, it has a 
potential energy of the form 

V{(j)) = ^m^(j)^ , (j)<Mp , (6.86) 

in terms of the canonically normalized field corresponding to its collective coordinate in 
the Ml direction. This collective coordinate will play the role of the inflaton, and we will 
refer to this D4-brane as the inflaton brane. 

As mentioned above, at mi = j/M,j = 1,...,M there is a rectangular torus in the 
M2, x' directions, equivalent by an SL(2,Z) transformation to the one at the origin. Introduce 
N4 additional D4-branes wrapped on the corresponding SL(2,Z) transforms of the cycle 
wrapped by the inflaton brane. The jth such brane has a quadratic potential proportional 
to (mi — j/My, minimized at Mi = j/M. Place each at its minimum. As the inflaton brane 



rolls down its potential (6.86), it encounters these additional branes, causing the strings Xj 



(and fermion partners) stretched between them to come down to zero mass. That is, (f) and 



the Xj couple as in our basic field theory model (1.1). 

It is clear that this structure arises more generally than the particular model [16,32]. In 
this particular case it is worthwhile to analyze the consistency of these added branes with 
the moduli stabilization barriers introduced by the curvature of the nilmanifold and other 
ingredients required to stabilize the space. A single D4-brane at the minimum of its potential 
is subdominant to the moduli-stabilizing barriers. There is a limit to how many additional 
branes can coexist with moduli stabilization. The tension of the set of D4-branes is 



V^^ = N. "Z^"" , . (6.87) 
(27r)4^,a'2 ^ ' 
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The directions ui and U2 are on the same footing; similar statements apply with the two interchanged 

M 
2 



and with x' replaced by x" = x + M-U1U2 
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This must be less than the scale of the moduli-stabilizing barriers, of order the curvature- 
induced potential energy: 

< ,^ " , . (6.88) 
Now in terms of the field theoretic quantities of the previous sections, A^4 ~ 0/A. So the 



condition (6.88) translates into the condition 



iV4 ~ < ——^ ^ . 6.89 

In the simplest version of the construction [32] - with the numerical examples discussed there 
and in [16] - the number of D4-branes is limited by this back reaction to be of order 10. 
Possibilities for warping down excessive contributions to the potential energy were discussed 
in [32]. In general, the mechanism we have discussed arises in a wide variety of "monodrofold" 
type compactifications [34]. 

A similar structure, with somewhat more flexibility in the parameters, arises in the 
setting [17] to construct trapped inflation from string theory. Consider type IIB string theory 
on a warped Calabi-Yau manifold, with an axion c arising from a 2-form RR potential C-^-* 
integrated over a 2-cycle S2. In the presence of an NS5-brane wrapped on S2 within a 
warped region (with a corresponding anti-brane wrapped on a homologous cycle in a distant 
warped region), the potential for c takes the form 



where e encodes the warp-factor dependence. As explained in [17], the axion decay constant 
is of order / ~ Mp/L^. This setup, with a large stabilized 2-cycle size naturally realizes 
large-field inflation (with c playing the role of the inflaton, executing many cycles of its 
basic period c — c + (27r)^). For the case of a blown-down 2-cycle, £ — 0, the same setup 
leads to trapping as follows. When £ = 0, as c rolls through the values (27r)^j (with j an 
integer), new light degrees of freedom appear in the theory. One intuitive way to see this is 
via the S- and T- dual setup depicted in Fig. 1 of [17] - there the NS5-branes' horizontal 
separation corresponds to £, and when this vanishes the unwinding motion takes the system 
through configurations where these NS5-branes meet. At these points, new light degrees of 
freedom arise from stretched D2-branes; the theory at low energies is a nontrivial interacting 
CFT (see e.g. [35]). The massless degrees of freedom of this CFT are produced much in 
the same way as are the x's described above (though perhaps in this case we should call 
it wnparticle production, since the low-lying degrees of freedom of the CFT are not strictly 
speaking particle states). In the original duality frame, the light "tensionless string" degrees 
of freedom arise with £ — from wrapped D3-branes (with appropriate worldvolume flux to 
cancel the contribution of cQs to the brane tension at the quantized values Cj = (2'7r)^j). 
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Of order 10^ to 10^ circuits can fit inside the compactification, satisfying tlie back reaction 
constraint delineated in eqn (3.42) of [17] by using the freedom to obtain somewhat large 
volume while maintaining high moduli stabilizing barriers using for example the methods of 
the large volume scenario [36] as explained in §4.4.1 of [17]. 

This construction corresponds to a linear potential V{(j)) = fi^(j) (modulated by instanton- 
generated sinusoidal corrections), a simple generalization of the m^0^ model analyzed above 
in §4. This potential is slightly flatter but leads to similar conditions on its parameters, 
shown in figure [3j 

Now for a potential V{(j)) = from our solution above we have 

2«B, , , C MpV""^-''"' / MA'"''-'" , , 



For m^0^, i.e. a = 2, this becomes (using (4.73)) 



A (27r)V/^iVe^^ V 

The number of e-foldings is 

]^Q-6 /^\ 4/3~a/3 



which reproduces the result (4.73) listed above for the case a = 2. 

The non-Gaussianity constraint m < lOif corresponds to (for a < 4) 



6 



For a = 2 this becomes 



and for a = 1 it is 



(6.92) 



< 1 . (6.94) 



m 10"-^° 

^ > {2.)-^g-"'lQ-'N:'l^ . (6.96) 
This corresponds to a constraint on the number of production events 

Nevents ~ ^ < {2T,f'^lQ~^N,g-^ (6.97) 

for any a. Similarly, one can derive a lower bound on Nevents from our slow roll conditions. 

For a wide range of parameters, the most stringent condition comes from (2.21). For a < 4, 
the constraint is 

Nevents > {27f)^/Ho-^^/^ N^g''' . (6.98) 
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Therefore, for any a < A one has a window of ~ 10^/^ between the minimum and maximum 
values allowed. 

Altogether, we find that the structure required for trapped inflation arises in the directions 
with monodromy in string compactifications, within a different regime of the potential and 
field range from that considered in modeling chaotic inflation in [16,17]. The ingredients 
required for trapped inflation generally introduce more back reaction than occurs in the 
corresponding single-field chaotic inflation model, but do fit into a reasonable subset of the 
known constructions. 

7 Discussion 

One of the satisfying recent developments in inflationary theory has been a more systematic 
classification of inflationary mechanisms. An inflationary mechanism can be characterized 
by its number of degrees of freedom - single field versus multiple field (a feature correlated 
with /nl^^ and isocurvature effects), the sound speed of its perturbations (correlated with 
/nl"^^^*'^'^'^^)' ^^"^ range of its inflaton (correlated with the gravity wave signature r). 

The present mechanism involves multiple fields (including the xs), but it behaves like some 
single field models in its prediction for large /^'^"'''^^'^'^^^j^ 

Although its signatures are somewhat similar to its strong-coupling analogue [6], we have 
seen that trapped inflation fits concretely into previously studied string compactifications; 
it is fair to say that the mechanism [5, 6] lacks a known clean top-down embedding (in 
the small subset of string compactifications yet studied). It would be interesting to find 
a compactification that interpolates between the two cases by varying the number of light 
degrees of freedom (and hence the 't Hooft coupling). 

The calculations in this paper required somewhat novel techniques for treating the effec- 
tive dynamics of resulting from the production of the sets of (temporarily) light x particles. 
There are several ways in which our analysis could be extended. In particular, it would be 
useful to develop more precise analytical tools to treat the perturbations. 
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